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Abstract. We present an experimental study of the flow dynamics of a lamellar phase sheared in the Couette 
geometry. High-frequency ultrasonic pulses at 36 MHz are used to measure time-resolved velocity profiles. 
Oscillations of the viscosity occur in the vicinity of a shear-induced transition between a high-viscosity 
disordered fluid and a low- viscosity ordered fluid. The phase coexistence shows up as shear bands on the 
velocity profiles. We show that the dynamics of the rheological data result from two different processes: (i) 
fluctuations of slip velocities at the two walls and (ii) flow dynamics in the bulk of the lamellar phase. The 
bulk dynamics are shown to be related to the displacement of the interface between the two differently 
sheared regions in the gap of the Couette cell. Two different dynamical regimes are investigated under 
applied shear stress: one of small amplitude oscillations of the viscosity {Si]/rj ~ 3 %) and one of large 
oscillations {Sr]/rj ~ 25 %). A phenomenological model is proposed that may account for the observed 
spatio-temporal dynamics. 
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1 Introduction 

When submitted to a shear flow, complex fluids exhibit 
unusual behaviour: steady shear modifies their structure 
and induces new structures or textures that do not exist 
at rest ^ . In wormlike micelles for instance, shear may in- 
duce a nematic phase j2l3j . whereas in membrane phases, 
shear creates new textural organizations such as multil- 
amellar vesicles j4l5j . Those shear- induced effects may be 
summarized in a shear diagram that specifies the nature of 
the stationary states (texture or phase) observed for differ- 
ent shear rates (or shear stresses). These states are sepa- 
rated by dynamical transitions that correspond to a jump 
between two branches of steady states In some cases, 
the fact that the complex fluid is sheared, and thus main- 
tained out of equilibrium, leads to richer behaviours than 
simple stationary states. Indeed, transitions to oscillating 
states and to chaotic-like behaviours have been reported 
recently in various systems (7i8-9m10.1L12.13.14) . These 
transitions involve rich dynamics of the fluid microstruc- 
ture, which arc often coined "rheochaos" ^Hl, and differ 
from classical instabilities, such as the Taylor-Couette in- 
stability or elastic instabilities )17) . 

In this work, we study a particular complex fluid ex- 
hibiting such structural dynamics: a lamellar phase com- 
posed of sodium dodecyl sulfate (SDS), octanol, and brine. 
The shear diagram of this lyotropic system depends upon 
the temperature. For small applied shear rates 7 < 71, 
the membranes are rolled into multilamellar vesicles called 



"onions" , that form a disordered close-compact monodis- 
perse assembly ^H] . The characteristic size of the onions 
is a few microns. For large applied shear rates 7 > 72 > 71 
and at temperatures T ^ 30 °C, onions are hexagonally 
ordered into layers that lie in the velocity-vorticity plane 
normally to the shear gradient and slide on one another. 
This structural transition, called the "layering" transition, 
is shear-thinning and occurs between 71 and 72 at a given 
critical shear stress CTc QHI' The rheological flow curve 
(measured shear stress a vs. applied shear rate 7) thus 
shows very strong shear-thinning around ctc that leads to 
a quasi-horizontal plateau at Cc- Such a rheological be- 
haviour seems generic in complex fluids as soon as a shear- 
induced transition is involved |19l20l2ll^ and has been 
much studied in the framework of shear banding theories 

mm . 

Under imposed shear stress^ the system follows the 
same rheological flow curve. Stationary states are reached 
at low and high shear stresses below and above ctc- How- 
ever, in a narrow region in the vicinity of the critical 
shear stress, complex temporal rheological behaviours are 
recorded and no stationary state is ever reached |10lllj . 



These temporal behaviours correspond to structural chan- 
ges of one part or of the entire sample between a disordered 
and an ordered assembly of multilamellar vesicles. 

In order to obtain good reproducibility in the flow 
curve a vs. 7, careful procedures must be followed: the 
applied shear stress must be maintained during time inter- 
vals 5t longer than 1000 s and the increment between two 
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different imposed shear stresses 5a must be smaller than 
1 Pa. However, the obtained state still depends upon the 
value of 5a and 5t. Using a rather fast procedure, we have 
evidenced in the plane cr vs. 7 a succession of stationary 
states, noisy regions, regions of small amplitude oscilla- 
tions ((57/7 = (577/77 ~ 3 % 121]) and of large metastable 
oscillations {Sj/'j ~ 25 %) as the shear stress is increased. 
A more "quasistatic" procedure reveals chaotic-like ape- 
riodic signals and the large oscillations are no longer ob- 
served, which confirms their metastability. Using dynam- 
ical system theory, we have shown that the chaotic-like 
oscillations of the viscosity may not be simply described 
by a three-dimensional dynamical system, probably be- 
cause spatio-temporal effects play a crucial role [TT] . 

In order to go further into the understanding of such 
behaviours, local velocity profiles are required to probe the 
flow inside the gap of a Couette cell (concentric cyhnders). 
Using dynamic light scattering (DLS) experiments in the 
heterodyne geometry, we have shown that a simple shear 
banding scenario holds for temporally averaged velocity 
profiles along the layering transition under applied shear 
rate the low-viscosity ordered phase nucleates at the 
rotor and progressively fills the entire gap of the Couette 
cell as the applied shear rate is increased from 71 to 72. In 
this range of applied shear rates, the two phases coexist 
and the flow is inhomogeneous. The temporally averaged 
position of the interface between the two phases lies at 
a fixed shear stress a*. However, noticeable fluctuations 
of the band position have been reported |2Z1- The poor 
temporal resolution of our DLS setup has prevented us 
to further study the dynamics under applied shear stress. 
Indeed, about three minutes are required to measure a 
complete proflle with the DLS technique |2H1- Note that 
this resolution is poor but comparable to that of Nuclear 
Magnetic Resonance (NMR) velocimetry with the same 
spatial resolution of 50 fim j29l3(J| . 

In this work, we present an experimental study of the 
layering transition under imposed shear stress using a re- 
cently developed experimental setup. High-frequency ul- 
trasonic pulses at 36 MHz are used to measure velocity 
proflles. Our technique is based on time-domain cross- 
correlation of high-frequency ultrasonic signals backscat- 
tered by the moving fluid. Post-processing of acoustic data 
allows us to record a velocity profile in 0.02-2 s with a spa- 
tial resolution of 40 fim |31| . Such a temporal resolution 
allows us to follow the dynamics of velocity profiles dur- 
ing the viscosity oscillations and to better understand the 
mechanisms at play during these oscillations. The article 
is organized as follows. In sectional we briefly describe 
the system under study as well as the experimental setup 
for measuring time-resolved velocity proflles in a Couette 
cell. Section 13 describes the experiments performed in the 
region of small amplitude oscillations. Scction^ldeals with 
the large metastable oscillations. Finally, Sec. [His devoted 
to a discussion of our experimental results and a phe- 
nomenological model is proposed that may account for 
the observed spatio-temporal dynamics. 



2 Experimental section 
2.1 System under study 

The complex fluid investigated in this work is made of 
sodium dodecyl sulfate (SDS), octanol, and brine. At the 
concentrations considered here (6.5 % wt. SDS, 7.8 % wt. 
octanol, and 85.7 % wt. brine at 20 g.L^^), a lamellar 
phase is observed at equilibrium. The smectic period is 
15 nm and the bilayer thickness is about 2 nm |32| . For 
the given range of concentrations, the lamellar phase is 
stabilized by undulating interactions. This system is very 
sensitive to temperature: for T > 35 °C, a sponge-lamellar 
phase mixture appears. 

Moreover, this lamellar phase is transparent to ultra- 
sonic waves and does not signiflcantly scatter 36 MHz 
pulses. Thus, in order to measure the velocity proflles, 
the sample has to be seeded with latex spheres of diam- 
eter of a few microns (see Sec. 12.21 below "l. These spheres 
are obtained from polymerization of an emulsion of di- 
vinylbenzene (DVB) and benzoyl peroxide (respectively 
39.6 % wt. and 0.4 % wt.) in water stabilized by SDS 
(1 % wt.) at a temperature of 80 °C 123- Benzoyl peroxide 
was added to the DVB to initiate the reaction. This emul- 
sion is then concentrated by centrifugation, diluted in the 
lamellar phase, and again concentrated. Finally, 1 % wt. 
of the obtained paste is suspended in the lamellar phase. 
Such a washing process allows us to avoid any modifl- 
cation of the lamellar composition. We checked that the 
addition of latex spheres does not signiflcantly affect the 
rheological behaviour and the layering transition of our 
system. 



2.2 Experimental setup 

To study the effect of shear flow on this lamellar phase, 
we used the experimental device sketched in Fig. ^a) and 
described at length in Ref. A rheometer (TA Instru- 
ments ARIOOON) allows us to perform rheological mea- 
surements in a Couette cell. The Couette cell is made of 
Plexiglas and has the following characteristics: inner ra- 
dius Ri = 24 mm, outer radius R2 = 25.07 mm, gap 
e ^ R2 — Ri = 1.07 mm, and height iJ = 30 mm. The 
walls of both cylinders are smooth. The rheometer im- 
poses a constant torque F on the axis of the inner cylin- 
der, which induces a constant stress a in the fluid, and 
measures the rotation speed f2 of the rotor (rotating inner 
cylinder) . A computer-controlled feedback loop on the ap- 
plied torque F can also be used to apply a constant shear 
rate 7 without any signiflant temporal fluctuations. The 
relations between the "engineering" {i.e. global) quantities 
((7,7) given by the rheometer and {F,f2) read 
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Such deflnitions ensure that (tr, 7) correspond to the spa- 
tially averaged values of the local stress and shear rate 
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in the case of a Newtonian fluid. The whole ceU is sur- 
rounded by water whose temperature is kept constant to 
within ±0.1 °C. 




Fig. 1. Experimental setup, (a) General sketch of the experi- 
ment, (b) Enlarged view of the gap of the Couette cell as seen 
from above. Dotted lines represent the ultrasonic beam and 
black dots (•) indicate the latex spheres that contribute to the 
backscattered signal. 



The thickness of the stator (fixed outer cyhnder) is 
2 mm everywhere except for a small rectangular window 
where the minimal thickness is 0.5 mm in order to avoid 
additional attenuation of the ultrasonic pulses due to Plex- 
iglas. A piezo-polymer transducer of central frequency / ~ 
36 MHz is immersed in the water in front of the stator. It 
generates ultrasonic pulses focused in the gap of the cell 
with a given incidence angle 9i ~ 10° relative to the nor- 
mal to the window. The ultrasonic pulses travel through 
Plexiglas and enter the gap with an angle 9 that is given by 
the law of refraction (see Fig.^b)). The transducer is con- 
trolled by a pulser-receiver unit (Panametrics 5900PR). 
The pulse repetition frequency is tunable from to 20 kHz. 
Typically, pulses are separated by a few milliseconds. Once 
inside the fluid, ultrasonic pulses get scattered by the latex 
spheres. Backscattered signals are sampled at 500 MHz, 
stored on a high-speed PCI digitizer (Acqiris DP235), and 
transferred to the host computer for post-processing. Typ- 
ical recorded backscattered signals are 1000 point long, 
which corresponds to a transit time of 2 /xs. 



Under the assumption of single scattering, the signal 
received at time t can be interpreted as interferences com- 
ing from scattering particles located at position y = cot/2, 
where cq is the sound speed in the fluid and y the distance 
from the transducer along the ultrasonic beam. When the 
fluid is submitted to a shear flow, the backscattered signals 
change as the scattering particles move along. Two suc- 
cessive pulses separated by AT lead to two backscattered 
signals that are shifted in time. The time-shift 5t between 
two echoes received at time t corresponds to the displace- 
ment 5y = co^i/2 along the y-axis of scattering particles 
initially located at position y = CQt/2. The velocity Wj,, 
projected along the j/-axis, of the scatttering particles lo- 
cated at y is then simply given by Vy = Sy/AT. A cross- 
correlation algorithm is used to estimate the time-shift St 
as a function of t. The sound speed cq is measured inde- 
pendently, which yields the projection Vy{y) in the entire 
gap. 

Finally, a calibration procedure using a Newtonian fluid 
is required to flnd precisely the angle 9, so that veloc- 
ity proflles v{x) of the tangential velocity v as a. func- 
tion of the radial distance x to the rotor can be computed 
from Vy{y) under the assumtion that the flow remains two- 
dimensional and axisymmetric. The reader is referred to 
Ref. for more technical details. Typical velocity pro- 
files are obtained by averaging over 50 series of twenty 
consecutive pulses. The spatial resolution of our experi- 
mental setup is of the order of 40 fxui and the temporal 
resolution ranges between 0.02 and 2 s per profile. 

3 Small amplitude oscillations 

3.1 Experimental procedure and rheological 
measurements 

To understand how the system flows when oscillations 
of the viscosity occur, we record the "engineering" (i.e. 
global) rheological data simultaneously to the velocity pro- 
files. In order to obtain reproducible experiments, care- 
ful rheological protocols must be followed as discussed in 
Ref. m]. At the temperature under study (T = 32 °C), 
we apply a first step of imposed shear stress at = 15 Pa 
during 6300 s. This step of applied shear stress allows 
us to start the experiment with a well-defined stationary 
state of disordered onions. We then apply increasing shear 
stresses for at least 2800 s per step. The stress increment 
between two steps is Sa = I Pa. 

Engineering rheological measurements are displayed in 
Fig. 121 In the imposed shear stress mode, the tempo- 
ral fluctuations of the shear stress are completely negli- 
gible. For the first two steps, the system reaches a sta- 
tionary state ((57/7 ^ 0.7 %). However, for imposed shear 
stresses greater than 16 Pa, noticeable fluctuations of the 
measured shear rate (and thus of the global viscosity of 
the sample) are recorded: S'j/j ~ 3 %. The character- 
istic times involved in such temporal behaviours are in 
the range 50-500 s as reported in previous experiments 
[101111 . Moreover, using static light scattering, we checked 
that the onset of these temporal fluctuations corresponds 



4 



S. Manneville et at: Spatio-temporal study of rheo-oscillations using ultrasound 



CO 



(a) 



40 



30 



« 20 



10 



(b) 














- 


1 ^ 









10000 20000 
t(s) 



Fig. 2. (a) Imposed engineering shear stress o-{t). Thick lines 
correspond to the time intervals when velocity profiles were 
recorded, (b) Measured engineering shear rate 7(t). The dotted 
lines correspond to the end of the various steps. The tempera- 
ture isr = 32 °C. 



to the apparition of the ordered texture in the sample 
(see Refs. |1UI11| for more details). Finally, note that all 
the experiments discussed here and in previous papers arc 
conducted far from any hydrodynamic or elastic instability 
(such instabilities would, in any case, lead to much shorter 
time scales) , so that the observed temporal behaviours can 
unambiguously be attributed to the presence of a shear- 
induced structural instability in our lamellar phase. 



3.2 Time-averaged velocity profiles 

In Fig. 13 we present a time average of all the tangential 
velocity profiles v{x) measured for the different applied 
shear stresses, x denotes the radial distance to the ro- 
tor (a; = at the rotor and a; = e at the stator). For 
each applied shear stress, the time intervals when veloc- 
ity profiles were measured are indicated by thick lines in 
Fig. Pl^a). Typically one hundred profiles (recorded every 
5 s in about 1 s per profile) have been averaged to obtain 
the data of Fig. 13 The error bars correspond to the stan- 
dard deviation of the velocity measurements. They are 
much larger than the uncertainty on individual measure- 
ments and are due to strong temporal fluctuations of the 
velocity field. These profiles clearly reveal an inhomoge- 
neous flow. Two differently sheared bands with different 
viscosities coexist in the gap of the Couette cell. The dot- 
ted lines in Fig. 13 indicate the average position of the 
interface between the two regions of different viscosities. 
When the stress is increased, the highly sheared band ex- 
pands across the gap. This picture is in agreement with 
the classical shear handing scenario |23I24| and with pre- 



vious measurements obtained on the same lamellar phase 
(without seeding particles) using DLS |26| . 




Fig. 3. Time-averaged velocity profiles for a = VJ (•), 18 

(□), and 19 Pa (o) recorded simultaneously to the rheological 
data displayed in Fig.|5] Error bars correspond to the standard 
deviation of the local velocities and mainly account for tem- 
poral fluctuations of the velocity field (see text). The dotted 
lines show the average position of the interface between the 
two shear bands. The mean velocity of the rotor is vq — 20.7, 
28.8, and 38.7 mm.s"^ respectively. Wall slip occurs both at 
the rotor and at the stator (see text). 



Moreover, significant wall slip can be detected on these 
measurements: the velocity at the stator does not vanish 
at a; = e and the velocity does not reach the rotor velocity 
vq = Rif2 at X = 0. Let us define the slip velocity at the 
rotor Vsi = vq — vi as the difference between the rotor 
velocity vq and the velocity vi of the fluid close to the 
rotor, and the slip velocity at the stator Vs2 = ^'2 as the 
velocity V2 of the fluid near the stator. Experimentally, vi 
and V2 are estimated by linear fits of the velocity profiles 
over 100 close to the walls. At the rotor, one gets Vsi = 
4.5, 4.1, and 3.8 mm.s~^ for engineering shear stresses a = 
17, 18, and 19 Pa respectively, which correspond to local 
shear stresses at the rotor cti ~ 17.7, 18.8, and 19.8 Pa 
respectively. At the stator, Vs2 = 1-4, 1.7, and 1.9 mm.s"^ 
for local shear stresses at the stator tT2 = 16.3, 17.2, and 
18.2 Pa respectively. The local shear stresses have been 
calculated using the following relationship: 



2i?? 



Rt 



(3) 



where j = 2 (resp. j = 1) when i = \ (resp. i = 2). Note 
that it is important not to confuse the engineering value 
a with the local values ui, <T2, and more generally '7{x) = 
aiRf/{Ri+x)'^. In the coexistence domain of the layering 
transition, we thus note a large difference between slip 
velocities at the rotor and at the stator: wall slip is much 
stronger at the rotor than at the stator. These results are 
in qualitative agreement with previous DLS experiments 
showing the existence of lubricating layers of thickness of 
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about 50-250 nm with larger slip velocities at the rotor 
than at the stator |2f)j . 

In the following section, we focus on the individual 
profiles measured at cr = 19 Pa and try to capture the 
details of the flow dynamics. 



3.3 Study of the step at ct = 19 Pa 

In this part of the article, we present a complete study of 
the last step at an imposed shear stress of 19 Pa. We chose 
to study this particular step because of the large number 
of velocity profiles recorded during this step (415 profiles 
separated by 5 s). Note, however, that the behaviour ob- 
served at 19 Pa is qualitatively the same for the two other 
steps where similar complex dynamics oi^{t) are recorded 
(cr = 17 and 18 Pa). 



3.3.1 Local shear rate and interface position 

Figure E^a) presents a spatio-temporal diagram of the 
shear rate profiles 7(0;, t) Pl]. The local shear rate is given 

by 

iM = -iR,^.)lJ^4^). (4) 



dx \Ri + X J 

Such a spatial derivative of the velocity field is very sensi- 
tive to experimental uncertainties and we used a moving 
average on a space window of three consecutive points and 
on a temporal window of three consecutive profiles to ob- 
tain the data shown in Fig.QJa). The abscissae in Fig.EIa) 
correspond to time t whereas the ordinates correspond to 
positions x in the gap. The value of the shear rate at a 
given time and a given position in the gap is coded using 
gray levels, darker values corresponding to smaller shear 
rates. Figure QJa) clearly reveals that the flow is both in- 
homogeneous and non- stationary. A highly sheared band 
is located near the rotor {x ~ 0-0.4 mm). The interface 
between the two shear bands is rather sharp and its posi- 
tion seems to fluctuate strongly in time. 

In order to present the evolution of the interface po- 
sition in more details, we calculate the derivative of the 
shear rate profile as a function of space d'y{x,t)/dx. In 
the case of a two-band velocity profile, one expects such 
a derivative of the local shear rate to be close to zero in 
the bands and to present a maximum in the region of the 
interface. Again, dj{x,t)/dx had to be filtered by using 
a spatial moving average over four consecutive points in 
order to reduce noise due to differentiation. The resulting 
data shown in Fig. ^b) reveal that the position of the 
interface between the two bands oscillates inside the gap. 
We also note that, in some particular cases indicated by 
white lines and dots in Fig. 2fb), two (or more) maxima 
are found at a given time. These two maxima correspond 
to two interfaces and thus to the presence of three bands 
of different shear rates. Such profiles with three (or more) 
shear bands have a short lifetime and seem highly un- 
stable: the interface near the stator either moves rapidly 
towards the interface located near the rotor or disappears 
at the stator (see below). 
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Fig. 4. (a) Local shear rate 7(3;, t). (b) Spatial derivative of the 
local shear rate : d-y{x,t)/dx. Moving averages on both space 
and time were used to filter out noise due to differentiation 
of the experimental data. White lines and dots indicate time 
intervals when two (or more) interfaces are present (see text). 



3.3.2 Individual two-banded and multi-banded velocity 
profiles 



Oo 
30 Pa» o 



I 20 
10 



□ • o 

□ • "o 




0.5 

X (mm) 



1 



0.5 1 
X (mm) 



Fig. 5. Individual velocity profiles recorded at cr = 19 Pa. (a) 
t = 617 (•), 863 (□), and 1000 s (o). The dotted lines show 
the position of the interface between the two shear bands, (b) 
t — 525 (□) and 536 s (•). The dotted line marks the position of 
the main stable interface whereas dashed lines help to locate an 
unstable secondary interface near the stator. The experimental 
uncertainty is of the order of the marker size. 



In order to better illustrate the two points discussed 
above (fluctuations of the interface position and presence 
of more than two shear bands in the Couette cell), we 
isolate different individual velocity profiles in Fig. [S| Fig- 
ure |3^a) unambiguously shows that the interface between 
the two shear bands fluctuates: it is located at x = 0.42 mm 
at < = 617 s, moves to x = 0.28 mm at t = 863 s, and 
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goes as far as x = 0.77 mm at t = 1000 s. These pro- 
files present a single interface whereas three interfaces are 
clearly seen on the velocity profiles shown in Fig.[5{b). At 
times t = 520-540 s, a very narrow, highly sheared band 
is nucleated near the stator and moves rapidly towards 
the rotor where it coalesces with the initial, more stable 
high-shear band located at a; ~ 0.4 mm. 

From the data of Fig.|SJb), we may also roughly esti- 
mate the velocity of the secondary interface to be about 
0.1 mm/10 s~ 0.01 mm.s~^. On the other hand, the 
main interface travels about 1 mm in typically 100 s (see 
Fig. Efb)), which corresponds to an interface velocity of 
the same order of magnitude. Thus, if one assumes that 
the propagation speed is characteristic of the interface 
[851 . we may infer that the various bands are of the same 
nature, although secondary bands seem unstable. Note 
that this picture assumes that the flow remains two-dimensii 
and axisymmetric so that the ultrasonic measurements 
can be interpreted in terms of the tangential velocity v{x, t). 
However, we cannot exclude the possibility of three-dimen- 
sional events that could lead to transient bumps in the 
velocity profiles such as those of Fig. |Sfb). Simultaneous 
imaging of two components of the velocity field is required 
to draw definite conclusions on that specific issue and is 
left for future work. 



3.3.3 Wall slip and band dynamics 

As shown above in Sec. 13.21 wall slip occurs both at the 
rotor and at the stator and the slip velocities Vsi and 
may well vary strongly in time. Thus, in order to analyze 
further the flow dynamics, we first need to measure quan- 
titatively the slip velocities as a function of time, as well 
as the shear rate in the bulk of the lamellar phase. 

Figure|SJa) presents the temporal evolution of Vsi and 
Vs2- In order to avoid any side effects, Vsi{t) and Ws2(i) 
are obtained from linear fits of each velocity profile for 
< X < 100 fim and e — 150 /Ltm< a; < e — 50 fim re- 
spectively. As noticed above on the average profiles, slid- 
ing is stronger at the rotor than at the stator. Like v{x, t), 
the two slip velocities present important temporal fluctua- 
tions with characteristic times in the range 50-500 s. Slip 
dynamics in the two lubricating layers are also strongly 
dissymmetric since the slip velocity at the rotor fluctuates 
by 5vsi/vsi — 37 %, whereas dvs2/vs2 — 15 % only. 

From Vsi and Vs2, and knowing the measured engineer- 
ing shear rate 7, we can compute the "true" shear rate 
7true produced by the lamellar phase and defined consis- 
tently with Eq. ^ as 

. Rl+Ii'2 VQ- Vsi ~ %Vs2 

i?i(i?i-fi?2) e ' 

where vq — RiQ is computed from 7 using Eq. Q (see also 
the Appendix). We recall that the measured engineering 
shear rate 7 is the sum of the contributions of the two 
lubricating layers at the walls and of the lamellar phase. 
As seen on Fig.lSJb), the true shear rate fluctuates much 
more than the engineering shear rate: 57truc/7truc = 5.5 % 




500 1000 1500 2000 

t(s) 

Fig. 6. (a) Slip velocities at the rotor Vsi{t) (thin hue) and at 
the stator 1)32 (i) (thick line), (b) Engineering shear rate 7(4) 
(thick Une) and true shear rate 7truo(i) (thin line). These data 
were filtered by using a moving average on a temporal window 
of three consecutive profiles. 



whereas 157/7 = ^vq/vq = 2.5 % in this case. Moreover, 
the engineering shear rate and the true shear rate may 
be completely out of phase due to the dynamics of slip 
velocities: 7true sometimes decreases whereas 7 increases 
as well as the slip velocities. This is the case for instance 
around t = 400 s and t = 1200 s. Thus, wall slip and 
bulk dynamics are decoupled under imposed shear stress, 
leading to the observed discrepancy between the 7 and 
7truc signals. Note that, if the experiments were conducted 
under applied shear rate (i.e. at fixed vq), the dynamics of 
7truc would be automatically coupled to that of Vsi + Vs2 
through Eq. (0). 

Figure presents the position of the main interface 
as a function of time together with the evolution of the 
true shear rate. The position S(t) of the main interface is 
defined as the position where d'j{x,t)/dx reaches a max- 
imum as a function of x. The black lines and dots in 
Fig. IJb) indicate the time intervals when this detection 
method may not be valid due to the presence of two inter- 
faces in the gap of the Couette cell. Two major conclusions 
may be drawn from these data. 

(i) First, the position of the interface 5{t) seems rather 
well correlated to the true shear rate 7truc (t) when a single 
interface is present in the gap (at least much better cor- 
related than to j{t)). When 7true(0 increases, the band 
moves towards the stator and tries to fill the entire gap. 
When 7truc(i) decreases, the band moves back to the ro- 
tor. This is consistent with a dynamical picture of shear 
banding where fluctuations of the interface position S(t) 
and of 7truo(0 would be linked instantly: an increase of 
7truc(i) corresponds to an increase of the proportion of 
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Fig. 7. (a) True shear rate 7truo(i)- (b) Position of the interface 
5 it). Black lines and dots indicate time intervals when two 
interfaces are present (see text). 
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Fig. 8. True shear rate 7truc(t) (thick line) and local shear 
rates 71 (i) (thin bottom line) and 72(f) (thin top line) in the 
weakly and highly sheared band computed from Eqs. 12311 and 
12411 respectively. The dotted lines correspond to estimations of 
71 (f) and 72(f) by averaging the local shear rate ■y(x,t) in the 
two bands. Black lines and dots indicate time intervals when 
the analysis is not strictly valid due to the presence of more 
than two shear bands (see text). 



the highly sheared band and thus to an increase of 5{t) 
(see also SecEJ. 

(ii) Second, such a correlation between 7truc {t) and (5(f) 
is sometimes far from perfect, even when only one interface 
is present (see t = 1600-1800 s), suggesting that the local 
shear rates in the two bands also fluctuate in time. Indeed, 
if the two local shear rates were constant and equal to 71 
and 72 in the weakly and highly sheared bands respec- 
tively, the continuity of the velocity field would impose 



that S{t) and jtmoit) arc linearly linked by 

. N . e-d{t) . S(t) . . , S(t) 

7t™c i - 71 ^ +72 — = 71 + 72 - 71 — , 6 

e e e 

at least in the presence of a single interface. 

Figure IHI presents the evolution of the local shear rates 
7i(f) and 72(f) in the two bands inferred from Vsi{t), 
Vs2{t), S{t), and v{S{t),t) according to Eqs. and 
(see Appendix). Another estimation of 71(f) and 72 (Oi 
obtained by averaging the local shear rate 7(x,f) over 
(5-1-100 /im< X < e— 100 fim and 100 ^m< x < (5—100 yum 
respectively, is shown in dotted lines in Fig. |S| The two es- 
timations of the local shear rates in the bands are in good 
quantitative agreement. In both bands, the local shear rate 
strongly fluctuates by 571/71 ~ 23 % and (572/72 — 12 %. 
Once again, the time scales involved in these fluctuations 
range from 50 to 500 s. Such large fluctuations imply that 
Eq. (jSJ with constant 71 and 72 does not hold in our ex- 
periment. Finally, note that 71(f) and 72(0 are clearly 
correlated but seem anti-correlated with (5(f). Moreover, 
the relative fluctuations of the local shear rates are about 
twice as large in the weakly sheared band as in the highly 
shear band (see Sec. for more details). 

In this first experiment, wc have analyzed the flow dy- 
namics in a regime of small amplitude oscillations. One 
may also wonder what happens when large amplitude os- 
cillations are produced such as those observed in Ref. ^U] ■ 
what are the similarities and the differences between those 
two dynamical regimes? In the following section, we show 
how to induce easily large amplitude oscillations and an- 
alyze the dynamics involved in such conditions. 



4 Large amplitude oscillations 

4.1 Experimental procedure and rheological 
measurements 

In this section, we study the regime of the large amplitude 
oscillations. This regime is difficult to attain experimen- 
tally by increasing the applied shear stress step by step as 
in standard procedures: in most experiments, just after the 
region of small amplitude oscillations, the system jumps 
directly to the highly sheared branch and the regime of 
large oscillations is missed. This seems to be related to 
the very small size of the coexistence region in the orien- 
tation diagram. Indeed, the large oscillations arc confined 
in a mctastable region of applied shear stress which is less 
than 1 Pa wide [10111) . Moreover, due to filling conditions 
of the Couette cell that may vary from one experiment to 
the other, the exact value of the shear stress that has to 
be applied to obtain such oscillations is not reproducible. 

In order to find this region more easily, wc first im- 
pose a constant shear stress of 12 Pa for one hour, fol- 
lowed by a constant shear stress of 13 Pa for another 
hour at a temperature T = 28 °C. This allows us to 
prepare a stationary state of disordered onions. We then 
bring the system into the coexistence domain by applying 
steps of constant shear rate at 7 = 15, 20, and 25 s~^ 
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walls. Indeed, as shown by the solid line in Fig. EH the 
velocity profile may be well fitted using a power-law fluid 
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Fig. 9. (a) Engineering shear rate 7(t). (b) Engineering shear 
stress (j{t). The dotted lines correspond to the end of the var- 
ious steps. The temperature is T = 28 °C. 



during about 4000 s. This procedure is easy because the 
coexistence zone is very large when the applied param- 
eter is the shear rate rather than the shear stress. We 
then note the measured shear stress (ct ~ 14.4 Pa at 
7 = 25 s""'^), switch the rheometer from the shear-rate- 
controlled mode to the stress-controlled mode, and apply 
a constant stress (cr = 14.4 Pa) for about three hours. 
Note that the ARIOOON rheometer allows us to program 
a new step in the rheological procedure without stopping 
the experiment. This allows us to begin the experiment 
in the coexistence domain at imposed shear stress, which 
induces the large oscillations of the viscosity shown in 
Fig. ini^a): a low viscosity {i.e. a high shear rate) is ob- 
served when the onions are in the ordered state whereas a 
high viscosity {i.e. a small shear rate) is measured when 
the onions are in the disordered state. The period of these 
large oscillations is about 2000 s and their amplitude is 
5r]/r] = (57/7 = 25 %. Finally, a last step at a = 14.2 Pa 
is performed, that shows that large oscillations are still 
observed at a slightly lower imposed shear stress. Figure|51 
shows the full rheological data sets 7(t) and a{t) recorded 
during such a procedure. 



4.2 Time-averaged velocity profiles 

Simultaneously to the rheological measurements, we record 
the velocity profiles. Figure ^| presents averaged velocity 
profiles obtained during pre-shearing at low imposed shear 
stress (ct = 12 Pa) and under imposed shear rate. When 
cr = 12 Pa, and thus at very low shear rate, the flow is 
homogeneous aird stationary. Moreover, the velocity pro- 
file is curved, revealing that disordered onions are highly 
shear-thinning. We also note that the onions slip at both 



= AY,, 



with 



0.2, as long as slip velocities are 



accounted for. 




Fig. 10. Time-averaged velocity profiles for cr = 12 Pa (□), 
7 = 15 (•), 20 (o), and 25 s~^ (■). The corresponding rheo- 
logical data are shown in Fig. |5] except for cr = 12 Pa. Error 
bars correspond to the standard deviation of the local velocities 
and mainly account for temporal fluctuations of the velocity 
field. The solid line corresponds to the homogeneous flow of a 
power-law fluid (shear-thinning exponent n = 0.2 and parame- 
ter A = 9.3) with wall slip (see text). The dotted lines show the 
average position of the interface between the two shear bands. 
The mean velocity of the rotor is i;o = 5.9, 14.7, 19.6, and 
24.5 mm.s~^ respectively. The horizontal solid line indicates 
vq — 5.9 mm.s^ . 



When the shear rate is applied and increased above 
7i, the flow becomes inhomogeneous. A high-shear band 
nucleates at the rotor and coexists with a low-shear band. 
Moreover, the velocity profiles as well as the position of the 
interface fluctuate, as shown by the large errors bars on 
Fig. EH localized near the mean position of the interface. 
These fluctuations induce small variations of the measured 
shear stress {6a j a ~ 3 %) as seen in Fig. |3Jb). The same 
behaviour is encountered for all the steps at imposed shear 
rate in the coexistence region between disordered and or- 
dered onions |27j . Increasing the shear rate in this region 
only changes the proportion of the two structures. More- 
over, the evolution of the average slip velocities as 7 is 
increased is as follows: Vs\ = 5.6, 6.5, and 6.5 mm.s~^ and 
Vs2 = 2.4, 2.0, and 1.4 mm.s'^ for 7 = 15, 20, and 25 s'^ 
respectively. 

Finally, once the rheometer is switched to the stress- 
controlled mode, the fluctuations of the velocity field get 
so large that a global time-average of the velocity profiles 
becomes meaningless. In the following section, we pro- 
vide a detailed analysis of the last step at an imposed 
shear stress of 14.2 Pa. Qualitatively, the same behaviour 
is found when a = 14.4 Pa is applied. 
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4.3 Study of the step at ct = 14.2 Pa 
4.3.1 Velocity profiles v{x,t) 




30 



20 



These homogeneous profiles are stable during at least 
twenty minutes and then suddenly become unstable in a 
few minutes. The crossovers between the two states are 
very fast compared to the period of time spent in the two 
homogeneous states. As in Sec. 13.3.31 we need to quan- 
tify the dynamics at the walls before trying to describe 
precisely the flow field during the crossovers. 



4.3.2 Wall slip dynamics 



2000 3000 
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Fig. 11. Spatio-temporal diagram of the local velocity v{x,t). 
A moving average over three consecutive profiles was used to 
smooth the experimental data. 



Figure lTTI presents the spatio-temporal diagram oiv{x, t) 
corresponding to 760 profiles recorded every 6 s at cr = 
14.2 Pa. These measurements reveal that the velocity pro- 
files oscillate between two homogeneous states: a high and 
a low viscosity state |S3 ■ Figure El shows two velocity 
profiles averaged in the high and low viscosity states. The 
time intervals on which the average was taken are indi- 
cated by horizontal bars in Fig. I13r b). In both cases, the 
profiles are homogeneous. Moreover, the lamellar phase is 
highly shear-thinning since the velocity profiles may be 
well accounted for by power-law fluids {a = Aj'^^.^^) with 
exponent n = 0.1. We can also note that the onions still 
slip at the two walls of the Couette cell, so that the average 
slip velocities had to be included in the fits of Fig. [T^ 




Fig. 12. Velocity profiles averaged over t — 1122-1577 s in the 
low-viscosity ordered state (o) and over t = 2172-2310 s in the 
high-viscosity disordered state (•) (see the horizontal bars in 
Fig. I13f b'l'). The solid lines correspond to homogeneous fiows 
of power- law fiuids (shear-thinning exponent ti = 0.1 in both 
cases and parameters A — 10.15 and 11.25 respectively) with 
wall slip (see text). 




Fig. 13. (a) Slip velocities at the rotor Vsi{t) (thin line) and 
at the stator Vs2{t) (thick line), (b) Engineering shear rate 
'y{t) (solid line) and true shear rate 7tiue(i) (dotted line). The 
horizontal bars indicate the time intervals on which the velocity 
profiles of Fig. 1121 were averaged. The parts in thick lines refer 
to the increase and the decrease of the shear rate described in 
Figs. 1141 and 1151 respectively. 



Figure El presents the slip velocities at the rotor and 
at the stator. Again, Vsi and Vs2 are significantly different. 
Fast fluctuations on time scales shorter than one minute 
are observed. Slower variations occuring on a time scale 
of about 1000 s are superimposed to the fast fluctuations. 
This long time scale is clearly correlated to the state of 
the lamellar phase: both slip velocities are larger (resp. 
smaller) when the lamellar phase is in the disordered (resp. 
ordered) state. 

Finally, using w<,i(i), Vs2{t), jit), and Eq. lO, we have 
calculated the true shear rate jtmeit) produced by the 
lamellar phase (see Fig. I13f b)'). Contrary to what hap- 
pens in the regime of small amplitude oscillations, 7truc 
oscillates in phase with 7 and its variations are very large 
compared to the variations of slip velocities. 
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4.3.3 Crossover from the weakly sheared state to the highly 
sheared state 
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Fig. 14. Crossover during the increase of the shear rate indi- 
cated by a thick line in Fig. ll3f b'l. (a) Local shear rate ^{x,t). 
White lines indicate the time intervals when the velocity pro- 
files presented in (b) were averaged. A linear gray scale is used 
to code the values of 7: black corresponds to 7 = 5 s~^ and 
white corresponds to 7 = 45 s^^. (b) Velocity profiles averaged 
over t = 251-298 s (o), 348-403 s (a), 465-526 s (□), 588-658 s 
(V), and 714-770 s (•). 



Figure IT^ shows that the system goes from the weakly 
sheared state to the highly sheared state via an inhomo- 
geneous flow. A high-shear band is nucleated at the rotor 
and progressively fills the entire gap. Indeed, when aver- 
aging the velocity profiles on about 60 s (as shown by the 
white lines in Fig.^Ja)), one can see that the local shear 
rates in the low-shear and high-shear bands first remain 
constant as the interface between the two bands expands 
towards the stator. This occurs for the three first profiles 
in Fig. ll4r b). which arc consistent with "frozen" snapshots 
of a classical shear-banded fiow at different increasing im- 
posed shear rates. 

However, the nucleation and displacement of the band 
does not occur smoothly during all the crossover. At t ~ 
500 s, the interface stops and the viscosity in the highly 
sheared band decreases whereas the viscosity in the weakly 
sheared band remains constant. Then, the interface moves 
towards the stator again and the viscosity in the high- 
shear band does not change anymore. At the end of this 
process, we are left with a homogeneous velocity profile. 

Even if such a process is not exactly reproduced from 
one crossover to the other, the general features of the tran- 
sition from the weakly sheared state to the highly sheared 
state remain the same for all the investigated crossovers: 
steady growth of the band at fixed 71 and 72 seems to 
alternate with periods where the interface remains fixed 
while the local viscosities vary. 



4.3.4 Crossover from the highly sheared state to the weakly 
sheared state 

However, the above picture does not hold anymore dur- 
ing the crossover from the low viscosity state to the high 
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Fig. 15. Crossover during the decrease of the shear rate indi- 
cated by a thick line in Fig. I13f b). (a) Local shear rate 7(3;, t). 
White lines indicate the time intervals when the velocity pro- 
files presented in (b) were averaged. A linear gray scale is used 
to code the values of 7; black corresponds to 7 = 5 s~^ and 
white corresponds to 7 = 45 s~^. (b) Velocity profiles averaged 
over t = 1686-1754 s (o), 1807-1855 s (a), 1903-1955 s (□), 
2003-2065 s (V), and 2139-2198 s (•). 



viscosity state. In this case, as seen on Fig. the fiow 
remains homogeneous across the Couctte cell as the shear 
rate decreases. Thus, onions become disordered simulta- 
neously everywhere in the gap. The viscosity progressively 
increases without any significant inhomogencities and the 
true shear rate decreases accordingly. 

Once again, even if inhomogeneous flows are sometimes 
observed transiently on individual profiles during the re- 
turn to the disordered state, this general picture holds for 
the various investigated crossovers from the highly sheared 
state to the weakly sheared state. 



5 Discussion and model 
5.1 Summary 

Let us now summarize the results obtained in the present 
work. Under applied shear stress and inside the coexis- 
tence region of the layering transition, the measured en- 
gineering shear rate oscillates and the flow is inhomoge- 
neous. Two bands (or even more) of different viscosities 
coexist in the gap of the Couette cell and their relative 
amount varies as a function of time. Wall slip is also very 
important. The recorded signal 7(<) is the sum of three 
separated dynamics: the dynamics of the two slip veloci- 
ties and the dynamics of the lamellar phase. 

In the bulk, the true shear rate 7true(0 varies mainly 
because the interface between the two differently sheared 
bands oscillates. To a lower extent, the variations of 7true(0 
are also linked to the fluctuations of the local viscosities 
{i.e. of the local shear rates) in the two bands. In some 
particular cases, a highly sheared band is nucleated close 
to the stator. This secondary high-shear band is unstable 
and either dies at the stator or moves rapidly towards the 
rotor where it coalesces with the initial highly sheared re- 
gion. Most of the above features are shared by both the 
small amplitude oscillations and the large amplitude os- 
cillations. The main difference between large oscillations 
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and small oscillations lies in the fact that large oscilla- 
tions require the displacement of the band across the en- 
tire Couette cell whereas during small oscillations, the flow 
is always inhomogeneous and the interface between bands 
fluctuates only in one part of the cell. Moreover, in the 
regime of large oscillations, a dissymmetry is observed be- 
tween the crossovers from one homogeneous state to the 
other: the apparition of the ordered state seems to follow 
a nucleation and growth scenario whereas its destruction 
occurs homogeneously across the gap. 



5.2 Discussion 

Let us discuss our results in the existing experimental 
and theoretical framework. First note that these results 
are in agreement with our former analysis of the rhco- 
logical data using dynamical system theory ^I]. In this 
work, we were unable to prove that the dynamics of 7(f) 
corresponds to a three-dimensional deterministic chaotic 
system even though it has several distinctive features re- 
sembling those of a low-dimensional chaotic system. We 
proposed that spatial degrees of freedom were responsible 
for this failure. The present experimental results clearly 
show that spatial degrees of freedom are indeed involved: 
the recorded dynamics involves (at least) three different 
dynamics that seem uncoupled to each other. However, it 
remains to be explained why so many different regimes 
are observed, some of them strikingly close to temporal 
chaos, and others looking much more turbulent. 

Second, let us see whether our experimental data arc 
consistent with the main theoretical approaches dealing 
with complex fluids under shear flow. In the last decade, 
various microscopic or phenomenological approaches as 
well as numerical simulations have been used to describe 
inhomogeneous flows of sheared complex fluids |36I37I38I39| 
In those works, shear banding is captured by invoking mul- 
tiple flow branches in the constitutive relation a vs. 7, cor- 
responding to different microstructures. Under shear, the 
system separates into a steady state composed of shear 
bands, each lying on its own flow branch. The presence 
of spatial derivatives in the dynamical equations of mo- 
tion selects a unique stress cr* at which the interface be- 
tween two shear bands may be stable |37l40l4lll^ . This 
explains the experimental degeneracy in the shear stress 
at which coexistence occurs and assumes that the inter- 
face between the two bands is stable only at cr*. In the 
Couette geometry, the stress inhomogencity implies that 
the highly sheared band nucleates at the rotor and is sep- 
arated from the weakly sheared region by a single inter- 
face. It also leads to the following relationship between 
the shear stress at the rotor cti and the interface position 
5 inside the gap |26| : 




(7) 



In a previous work, we have shown that such a shear 
banding scenario holds for the layering transition in our 



lamellar phase [21] • It describes quantitatively the time- 
averaged profiles under both applied shear rate and shear 
stress (at least in the region of the small amplitude os- 
cillations). However, this oversimplified picture does not 
include the possibilty of oscillations of the rheological pa- 
rameters and of sustained oscillations of the position of 
the interface between the shear bands. Wc have proposed 
that these oscillations may be captured by assuming that 
the selected stress a* is a fluctuating parameter in our 
system Under the assumption of a fluctuating cr*, we 
may describe at least qualitatively the dynamics. Indeed, 
the same small fluctuations of a* will produce small vari- 
ations of the measured shear stress in the case of applied 
shear rate but huge variations of the measured shear rate 
in the case of applied shear stress. Thus, a microscopic 
equation was still missing to describe the evolution of a*. 

In the following, we first focus on the nonlocal ap- 
proach of shear-banded flows developed in Refs. )4()I4H . 
We then propose a simple set of equations based on such 
an approach but including the possibility of temporal dy- 
namics. This crude dynamical model may provide an ex- 
planation for variations of a* and for the observed dynam- 



5.3 The nonlocal approach of Refs. 



mm 



In Refs. |40l41j . it is assumed that the stress in the sys- 
tem is locally the sum of two terms, one related to the 
shear rate and the other one proportional to the second 
derivative of the shear rate: 



3(7) 



(8) 



The first term corresponds to the usual constitutive equa- 
tion for a homogeneous flow. The friction 17(7) between 
two adjacent sliding layers depends only on their relative 
velocity i.e. on the local shear rate. However, when large 
spatial variations of the local shear rate occur, a second 
term must be added that involves d^^/dx^. In this lat- 
ter case, the microstructure of the fluid varies on a length 
scale of the order of the range of interactions between 
the mesoscopic entities of the system (i.e. the microstruc- 
ture varies on the scale of a hydrodynamic layer). To take 
into account these inhomogenities of the microstructure, 
Dhont has introduced the "shear curvature viscosity" k 
(see Eq. ©) Ej. 

Like the standard shear viscosity, such a shear curva- 
ture viscosity depends upon the shear rate. In particu- 
lar, since the microstructure generally reaches a homoge- 
neous state at high shear rate, k has to decrease to zero 
when 7 goes to infinity. For example, for a fluid of rodlike 
molecules, shear progressively aligns the rods in the flow 
direction. At high shear rate, the rods are fully aligned 
and the microstructure cannot evolve further. Thus, the 
fluid becomes more homogeneous and the shear curvature 
viscosity has to vanish. To describe this behaviour, the 
authors propose that 



D 



1 -I- a 7^ 



(9) 
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where D is analogous to a diffusion coefficient and a is a 
phenomenological coefficient |4()I41| . 

Let us now analyze the consequences of such a be- 
haviour in the presence of a shear-induced transition. In 
this case, the first term 5(7) is a multivalued function of 
the shear stress: for a small range of applied shear stresses, 
two branches of different structures exist. In a plane geom- 
etry, Eqs. (jHJ and Q select a unique stress a* for which 
the two different structures coexist. If the applied shear 
stress is smaller than this selected stress a* , the system 
remains on the low-shear branch. Above a* , the struc- 
ture changes and the system jumps on the highly sheared 
branch. The value of a* is found by solving the following 
system of equations |4UI41| : 

r'K-ff(7))(l + «7')d7 = 0, (10) 

CT* = g{ii) = .9(72) , (11) 

where the unknowns 71 and 72 are the local shear rates in 
the low-shear and high-shear bands respectively. 



5.4 Toy model 

The nonlocal approach described above allows a good de- 
scription of shear-banded flows but is intrinsically sta- 
tionary and unable to reproduce any of the dynamics re- 
ported in the experiments. Recently, various phenomeno- 
logical models have been proposed to mimic the chaotic- 
like spatio-temporal behaviours of sheared complex fluids, 
mostly in the context of wormlike micelles and colloidal 
pastes |15I4HI44| . as well as liquid crystalline polymers 
|4RI48I47| . 

Here, we adapt the model of Refs. |40I41) by taking 
into account the particularity of our system, namely an 
assembly of soft multilamellar vesicles (onions) whose in- 
terlayer distance may vary under shear. Indeed, using neu- 
tron scattering, it has been measured that the smcctic pe- 
riod of the ordered onion texture decreases as the shear 
rate is increased whereas it remains constant in the dis- 
ordered texture It was assumed that some water is 
expelled from the onions and lies between the different 
planes of ordered onions to lubricate the layers. The time 
scale involved in such a release of water depends on the 
temperature, due to the presence of thermally activated 
defects in the lamellar phase. It ranges from ten minutes 
to an hour. Those water layers dramatically affect the vis- 
cosity of the sample. 

In the framework of Refs. |40l41j . such water layers 
should also strongly influence the coefficient a in the shear 
curvature viscosity (see Eq. lO). A scenario leading to 
sustained oscillatory behaviour then takes shape. Indeed, 
a decrease of a, linked to the slow release of water, may 
induce an increase in a* . Consequently, an increase of a* 
may lead to the destruction of the shear-induced phase, 
decrease the amount of water between the onions, and thus 
increase a (see below). The competition between those two 
processes may induce oscillations of a* and thus lead to 



oscillations of the measured shear rate under applied shear 
stress. 

More precisely, in order to mimic these two antago- 
nist effects, we introduce h{t) the total amount of water 
lubricating the planes of ordered onions. For a station- 
ary state in the coexistence region, we assume that h is 
directly proportional to the amount of ordered onions e. 
When the flow is composed of a single interface separat- 
ing a highly-sheared ordered region from a weakly-sheared 
disordered region, one simply has e = 5/e, where 6 is the 
position of the interface inside the gap of width e. Wc 
then write two coupled phenomenological equations obey- 
ing the following rules. 

(i) h must relax to /icq = in the absence of the shear- 
induced structure, to some maximum value /leq = when 
the fluid is fully ordered, and to the intermediate value 
/icq = e/io in the coexistence region. 

(ii) a must decrease when h increases. Indeed, when 
water is released, i.e. when h increases, the microstruc- 
ture of the ordered phase differs more sharply from that 
of the disordered phase. This should induce an increase of 
the shear curvature viscosity k, which in our model corre- 
sponds to a decrease of a. Thus, we assume that, at equi- 
librium, a relaxes to ao—f3heq, where ao and /? are positive 
parameters. Assuming exponential relaxations with char- 
acteristic times Ti and T2, wc propose to describe these 
two processes by 

§ = l(e(/)/i„-/i(/)), (12) 
at Ti 

^ = - (ao - /3/i(/) - a(t)) . (13) 
at T2 

Moreover, if Eq. {Tj) holds at all times, then e{t) is given 

where the stress at the rotor cti is linked to the applied 
shear stress by Eq. At each time step, cr*(/) is found by 
solving Eqs. (|10|I - H11|I using the current value a{t). The 
instantaneous relation H14|l between e{t) and cr*{t) sim- 
ply means that the time scales involved in the selection 
of CT* and in the transport of the shear stress are short 
compared to the swelling and deswclling processes of the 
onion phase. A more realistic approach would involve a 
third dynamical equation coupling 5 and cr* and/or mem- 
ory effects. This could lead to more complex dynamics 
than the simple oscillatory state found here, possibly to 
chaotic dynamics |15| . However, such an extension of the 
model is beyond the scope of the present article. 

In order to solve Eqs. (|10|1 - H14|I . a modelization of the 
rheological curve (7(7) describing the homogeneous part of 
the flow is required. For the sake of simplicity, wc assume 
that the two branches are Newtonian and we define 3(7) 
as 

5(7) = '717 if 7 < 7c, (15) 
.9(7) = mi if 7 > 7c, (16) 
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where 7c is some critical shear rate and 771 and 7/2 arc 
the viscosities of the disordered and ordered states respec- 
tively (771 > 772). Thus, for a G ['727c; ?7i7c], two different 
shear rates are accessible to the system. Eqs. ()15|l and 
H16|) provide a coarse approximation of the experimental 
flow curve but allow us to solve analytically Eqs. (|10|l - 
1)11(1 and to implement easily the numerical resolution of 
Eqs. inj-dl). 
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Fig. 16. Solutions of Eqs. (O-O for 771 = 0.4 Pa.s, 7)2 = 
0.1 Pa.s, 7c = 30 s~\ n = r2 = 200 s, ho = 30, ao = 0.2, and 
P = 0.024. (a) Selected shear stress (J*{t). (b) Global shear rate 
7(f). (c) Position of the interface 5{t). The imposed shear stress 
corresponds to a shear stress at the rotor ai = 5.2 (dashed 
line), 6.0 (solid line), and 6.5 Pa (dotted line). The dimensions 
of the Couette cell are _Ri = 24 mm and R2 = 25 mm. 



equivalently Eq. I|17|) ). and more generally, to assume that 
Eqs. (|5|l- (|ll|l are not significantly altered by curvature 
effects. 

Figure 1161 presents the temporal signals obtained us- 
ing the following set of parameters: 771 = 0.4 Pa.s, 772 = 
0.1 Pa.s, 7c = 30 s-i, ri = ra = 200 s, /iq = 30, ao = 0.2, 
and (3 = 0.024. Time series corresponding to (J*{t), 7(<), 
and d{t) are shown for various imposed shear stresses cor- 
responding to (Ti = 5.2, 6.0, and 6.5 Pa. For a small value 
of the applied shear stress (ai = 5.2 Pa, dashed line), the 
system relaxes towards a homogeneous stationary state. 
The system remains in the low-shear state since the shear 
rate is simply 7 ~ (Ti/771 ~ 13 s~^ and (5 = 0. For large 
applied shear stresses (ci = 6.5 Pa, dotted line), another 
homogeneous state is obtained but the viscosity of the sys- 
tem is now 772 since 7 ~ CTi/772 — 65 s~^. The fluid is then 
fully ordered (5 ~ e). 

For intermediate values of the applied shear stress (cJi = 
6.0 Pa, solid line), the two structures coexist and our sim- 
ple model predicts weakly damped oscillations. The pe- 
riod of the oscillations depends on the values of Ti, T2, 
P, and da* /da. The interface position 5{t) oscillates in 
phase with 'j{t) and in phase opposition with (7*{t) as well 
as 71 (t) = a* (t) 1 771 and 72 (i) = cr* (t) j 772 • Experimentally, 
the same behaviour has already been briefly mentioned in 
Sec. 13.3..'^ In Fig. El we plotted the experimental signals 
(5(t), 71 (i), and 72(0 normalized by their standard devia- 
tions and arbitrarily shifted vertically for better compar- 
ison. It is clearly seen that is indeed anti-correlated 
with 71 (i) and 72(^)7 which provides a nice confirmation 
for underlying temporal variations of a* . 
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Finally, the local shear rates 71 (t) and 72(0 are re- 
covered using Eq. Hll|l and the global shear rate '^(t) is 
computed by 



7truc(0 =7lW + (72(t) 



71 W) — 



(17) 



Note that, in principle, Eqs. l|S|l- l(TT)) and ((T7jl only hold 
in a plane geometry, whereas Eq. (|14|l results from the 
stress inhomogcneity inherent to the Couette geometry. 
However, as shown in the Appendix, the aspect ratio e/i?i 
of our Couette cell is small enough to use Eq. ^ (or 



Fig. 17. Normalized versions of the temporal signals of Figs. El 
and|Hl(see text), (a) Position of the interface <5(f). (b) Local 
shear rate in the weakly sheared band 71 (i). (c) Local shear 
rate in the highly sheared band 72(f). 



Note that, in the model, a*(i^ fluctuates by about 
1.2 % whereas large oscillations of 7(f) are found ((57/7 ~ 
25 %), which is consistent with the experimental data. 
However, due to the Newtonian assumptions of Eqs. I)15|l 
and the local shear rates 71 (t) and 72 (i) fluctu- 

ate by the same amount as G*(t) (~ 1.2 %), which con- 
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tradicts the experimental results of Sec. 13.3.31 that show 
much larger fluctuations. However, if nonlinear relation- 
ships 5(7) = Aij""' were considered, with i = 1 (resp. 
1 = 2) for the low-shear branch (resp. high-shear branch), 
Eq. ifTT)) would yield 6a*/a* = ruS^i/^i. This would lead 
to fluctuations of 7^ (i) enhanced by a factor l/rii as com- 
pared to those of <J*{t) and, if n\ ^ ni-, to some dissym- 
metry between the fluctuations of 71 (t) and 72 it) as seen 
in the experiments. Finally, note that the present model is 
also tractable under applied shear rate. Such an analysis 
is left for future work. 



Appendix: Determination of the true shear rate 
7truo aid of the local shear rates 71 and 72 

Let us consider the case of a Newtonian fluid of viscosity 77 
sheared in the Couette geometry. Moreover, let us assume 
that the fluid slips at both walls so that v(x = 0) = Ui 
and v(x = e) = V2- Using the local rheological behaviour 
cr{x) = Tj'^ix) and 



r 



2t:H{Ri+xY ' 



(18) 



6 Conclusions 

In this paper, we have shown that our experimental setup 
based on ultrasonic velocimetry allows us to follow the 
dynamics of the inhomogeneous flow of a complex fluid. 
We have defined precise rheological procedures in order to 
produce oscillations of either small or large amplitudes. 
Our experimental data unveil two major topics: 

(i) The dynamics observed on engineering rheological 
quantities is the sum of three a priori uncoupled systems: 
the two lubricating layers at the walls and the lamellar 
phase. This point is very important since it confirms the 
conclusions reached after an independent data analysis us- 
ing dynamical system theory |llj . In the present study, 
we have confirmed that the recorded signals cannot result 
from a low-dimensional deterministic chaotic system since 
(at least) three spatial systems are involved. 

(ii) The bulk dynamics is mainly due to the displace- 
ment of the interface between two differently sheared bands. 
Inspired by the ideas developed in Refs. |4t)l41| . we pro- 
posed that such a motion may be linked to the slow evo- 
lution of the selected stress a* . Assuming that the shear 
curvature viscosity depends upon the amount of water re- 
leased by the ordered onions, we have shown that a* could 
oscillate. These oscillations of a* lead to oscillations of the 
interface position 5{t) that are in phase opposition with 
the local shear rates 71 (t) and 72(0- The above features 
of the model are nicely confirmed by the experimental sig- 
nals, even though their dynamics are more complex. 

Note that the same kind of arguments may be invoked 
to describe the dynamics of the lubricating layers. More 
elaborate versions of the equations proposed here may in- 
duce complex behaviours closer to those recorded experi- 
mentally. 

Future experiments will focus on the measurement of 
the smectic period under imposed shear stress in the co- 
existence region, in order to reach a fully microscopic de- 
scription of the dynamics. Moreover, using cell walls of 
controlled roughness may help us to better understand 
the origin and the influence of wall slip dynamics. 

The authors wish to thank the "Cellule Instrumentation" at 
CRPP for designing and building the mechanical parts of the 
experimental setup. We are very grateful to A. Aradian, L. 
Becu, C. Gay, D. Roux, and A.-S. Wunenburger for fruitful 
discussions. 



where F is the torque applied on the rotor axis and H the 
height of the Couette cell, together with 



'^{x,t) = -{Ri+x) — 



d ( v{x,t) 



Ri 



one easily gets the following velocity profile: 



v{x) = V2 



Ri 



(Ri+x) 



r 



R2 "'AmjHRl 
Taking a; = in Eq. (IJOJ yields 



R2 



Ri 



Vi 



Ri 
R, 



r 



■ V2 



A-KTjHRiR^ 



{Rl 



Rl). 



(19) 

1 . 

(20) 

(21) 



Finally, with e = R2—R1, substituting a defined by Eq. 
into Eq. leads to 



Ttruc 



Rl 



Rl 



V, - ^^V2 



RiiRi + Ra) 



(22) 



which is exactly Eq. jSj) where v^i = vq —vi and Vs2 = V2- 
In the presence of a two-banded flow, we may now de- 
fine the local shear rate 71 in the weakly sheared band 
from the position 5 of the interface, the velocity v{6) mea- 
sured at X = S, and the slip velocity Vs2 at the stator as 



71 



{Rl + 5f 



Rl 



v{5) 



Ri+S 
R2 



Vs2 



{Ri+5){Ri+R2 + 5) 



(23) 



Equation (|23|l was simply obtained by taking i?i Ri+S, 
e — > e — 5, and vi — > v{d) in Eq. H22|l . Similarly, we define 
the shear rate 72 in the highly sheared band as 



72 



Rl + {Ri+S)^ VI - j^vjS) 
Rii2Ri + S) S 



(24) 



Combining Eqs. (j^ . (j^ . and (|^ . it is easy to see that 
Eq. ((BJ is not a strict equality but rather results from 
the small-gap approximation e/Ri <^ 1. Actually, using 
the dimensions of our Couette cell, we checked that the 
estimation of 7true given by Eq. 10) differs from that of 
Eq. (|22|1 by at most 4 %, which is of the same order of 
magnitude as the experimental uncertainty on the velocity 
measurements. 
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Finally, let us mention that our fluid is far from being 
Newtonian since shear-thinning exponents n = 0.1-0.2 are 
reported in Sec. 0] Repeating the above analysis with the 
local rheological behaviour cr{x) — ^7(2;)" leads to 



7tr 



2l-l/« ^Rl^Rlf/r. 



Ri 

Vl - — W2 
^2 



(25) 



instead of Eq. However, for our Couette cell, the 

difference between the estimation of Eq. 1)22(1 and that of 
Eq. H25|l is only 3 % in the worst case n — 0.1, which 
allows us to use Eqs. (|^ - lf^ to define 7truc, 7i, and 72, 
even if the fluid is not Newtonian. 
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